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2^ Abstract 

(N 

>■ We derive the two loop expressions for polygonal Wilson loops by starting from the one 

loop expressions and applying an operator product expansion. We do this for polygonal 
Wilson loops in R 1 ' 1 and find a result in agreement with previous computations. We also 
discuss the spectrum of excitations around flux tube that connects two null Wilson lines. 
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1 Introduction 

In this article we provide a short derivation for two loop polygon Wilson loop expectation 
value (or MHV amplitudes) using the Operator Product Expansion for Wilson loops in- 
troduced in [1]. This is done in a restricted kinematical region where the Wilson loop is 
embedded in an R 1 ' 1 subspace of four dimensional Minkowski space. 

The idea is very simple. One starts from the one loop result and then performs the OPE 
expansion. Each term of this expansion corresponds to the exchange of a free particle. Next 
we include the one loop energy or anomalous dimension for each particle. The expansion 
breaks the cyclic symmetry of the answer. If one completes it in the simplest way one 
obtains a two loop answer that has the correct expansion in all the OPE channels. The 
computation is particularly simple because at two loops one can organize the states into 
SL(2,R) representations. With R 1,1 kinematics only one SL(2) representation contributes 
(for each momentum), hence all states have the same anomalous dimension. 
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The two loop functions obtained in this fashion agree with explicit results computed 
by more direct methods [2] [3] or by other inspired guesswork [4]. These were checked 
numerically in [4] against the direct perturbative computation. 

This paper is organized as follows. We first review the OPE expansion for Wilson loops 
and explain why we have an SL(2) symmetry at two loops. We then discuss the derivation of 
the two loop results for the octagon and decagon from this points of view, leaving the general 
n case for appendix A. We then say a few words about the three loop octagon and present 
some conclusions. In appendix B we discuss some aspects of the spectrum of excitations 
around the flux tube that stretches between two null Wilson lines. This can also be viewed 
as the spectrum of excitations around the infinite spin limit of finite twist operators, or the 
GKP [5] string. 



2 OPE review 



The Operator Product Expansion for polygonal null Wilson loops is a certain expansion for 
the expectation value of a Wilson loop correlator in the limit where two or more consecutive 
edges become collinear. For more details see [1]. 

We consider a polygonal Wilson loop made out of n null segments. The collinear limit 
is approached as follows. We choose two null lines that belong to the Wilson loop polygon. 
These null lines are then extended. The two null lines preserve an SL{2) x R a x SO {2)^ 
subgroup of the conformal group. Using these null lines and two other arbitrary null lines 
we define a reference "square". This selects a dilatation operator R T inside SL(2), see figure 
1. This reference square is invariant under three commuting symmetries R T x R a x SO (2)$. 
Two of the symmetries are non-compact and are particular elements of the conformal group. 
Then a family of polygons is constructed by acting on the bottom part of the original polygon 
with these symmetries and joining it to the original top part of the polygon. Schematically 
we have 

(W)(T,a,<j)) = (top |e- TB+iffP+ ^ J (bottom) (1) 

This has an expansion in terms of intermediate states which are excitations of the flux tube 
that goes between the two selected null sides. Such states can also be described in terms of 
excitations around the infinite spin limit of local operators with spin. This correspondence 
with local operators enables the computation of the properties of the propagating states. 

A slight complication is that, strictly speaking, the Wilson loop correlation is zero. That 
is due to certain well understood UV divergencies. After introducing a UV cut-off one 
remains with a finite answer. These symmetries are broken in a precise way determined by 
an explicitly known anomaly [6]. We can thus take into account the effects of the anomaly. 

A standard way to take it into account is to consider the ratio function 

e« = (W/) r (2) 

[WmJ ' 

where the denominator is the same Wilson loop but in a U(l) theory, with the coupling 
replaced by T cusp . This ratio is a finite quantity, free of anomaly, known as the remainder 
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Figure V. We start from a general Wilson loop. We select two null lines, denoted in red. 
We then act with symmetries that are preserved by these null lines on one of the sides. The 
two null lines preserve a certain SL(2) symmetry. Selecting a dilatation operator inside this 
SL{2) amounts to a choice of a reference square, denoted here by the red dashed lines. This 
allows us to perform an expansion which has the rough form seen on the right. The first 
term comes from the exchange of the flux tube vacuum. The second from the exchange of 
a single excitation on the flux tube, the third from two excitations, etc. The denominators 
can be viewed as normalization factors, as we usually have in the standard OPE. 

function [6]. Unfortunately, in the OPE expansion of R one has already subtracted off some 
interesting contributions from the OPE expansion of {W)u(i), which is non-trivial. This 
makes it difficult to apply our program to directly to R. 

For us, in order to perform the OPE, it is convenient to introduce the following ratio 
function 



where the "square", "top" and "bottom" polygons are defined in figure 2. We can get rid of 
the divergencies completely if we choose a reference square that coincides with some of the 
cusps of the original polygon, see figure 2. 1 In this fashion we can have a well defined and 
non- vanishing answer already at one loop. The price we pay is however that r is not cyclic 
invariant. The remainder function differs from r in a simple way [1] 



and only the last term contributes to the OPE expansion. 

As an example consider the octagon in R 1 ' 1 kinematics, see figure 3. a. We label the 
location of the cusps as 



1 For a general reference square, the ratio (3) does not quite get rid of divergencies, it leaves some 
unimportant single logarigthmic divergencies which have zero momentum (under Here we will make 

the choice in figure 2 which gets rid of this ambiguity. 




(3) 



R — r — i?top + -^bottom ~~ ^cusp r U(l) ■ 



(4) 



{. . . , (x i , (Xj , x i ), (x i+1 , x i ),...} , 



(5) 




Figure 2: Definition of a ratio function that is finite and conformal invariant. It involves 
the selection of a reference square whose vertices (across a diagonal) coincide with vertices 
of the original polygon. 



where x — x ± t. For the Octagon, we fix these locations at 

xt = {-oo, -1, 0, e 2T } , xT = {-oo, -1, 0, e- 2CT } . (6) 
The corresponding ratio function at one loop reads 

r^° n = - y log (1 + e-^) log (l + e^) (7) 

This is the "seed" that we will use to get the two loop answer. This one loop result is the 
same as the result we would obtain in a £7(1) theory and it only comes from a single gluon 
(or single photon) exchange between the various null segments. 



2.1 SL(2,R) symmetry at two loops 

When we choose two null lines we preserve an SL(2) x R a x SO (2)^ symmetry which is 
a subgroup of the full conformal group, 5*0(2,4). Here we will make some remarks about 
the SL{2) symmetry. This SL{2) symmetry is broken by the flux tube. However, the mere 
existence of a flux tube is an effect of order g 2 . In fact, we will see that some consequences 
of the SL(2) symmetry are still preserved at low enough orders in perturbation theory. Let 
us start with the one loop answer. This arises from the propagation of free particles between 
the edges of a square, as in figure 1. A factor of g 2 already rises by creating and annihilating 
the excitation from the Wilson contour. Thus, the propagation in the bulk occurs as in the 
free theory. In the free theory this SL{2) is, of course, a symmetry, since it is just a subgroup 
of the full conformal group. Thus, the particles being exchanged form multiplets under this 
symmetry. An SL(2) representation consists of a primary P, or lowest weight state, and an 
infinite set of descendents of the form P n = (L_i) n P. More explicitly, in the case of R 1 ' 1 
kinematics, the only relevant operator is F + _ and the descendents are _D".F + _. The twist of 
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(a) 



(b) 




Figure 3: (a) The Octagon null Wilson loop in blue embedded in the Penrose diagram of 
R 1,1 . The red dashed line is the reference square we start with, (b) A different choice of 
reference square that is suitable for an OPE expansion in the same channel as in (a). In 
that example, the two choices of reference squares only differ by the position of the "bottom" 
cusp at (a) and — e (b). The corresponding OPE expansion parameters are related by the 
infinitesimal transformation r — > r + | e~ 2r . We must be able to re-write the OPE expansion 
after this transformation in the same form as before with the same anomalous dimensions. 
That is only possible if 7fc(jo) is independent of k. In (c) and (d) we draw the same picture 
as it would have looked like if the Octagon under consideration was Lorentzian. 



these operators is eo = 2 + 2n. The D + derivatives are essentially taken into account by the 
momentum quantum number, related to the R a . 

Now let us see what happens at two loops. At two loops, the particle starts feeling 
the effects of the color electric flux between the two null lines. This flux breaks the SL(2) 
symmetry. So, if we act with L_i on the flux vacuum, we can create a particle, but only with 
amplitude g, since the background flux F^ v ~ g 2 . In addition, the anomalous dimensions of 
the n-th single particle descendent, P n could be changed, to e = 2 + 2n + g 2r y n + o(g 4 ) . It 
is possible to show that 7„ is still independent of n, so that we have 

e = 2 + 2n + g 2 1 { P ) + o{g A ) , ra = 0,1,2,..- (8) 

where we noted that 7 can depend on p, which is the quantum number under the translation 
R a . This can be shown by starting with a single particle state P n \0). Acting with L_\ 
we get P n+ i|0) + gP n Po\0), where the second term is a two particle state arising from the 
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non-invariance of the vacuum. By acting with both sides of the commutator L_i = [L , L_i] 
on P n \0), we can show that 7 n should be independent of n. 

It turns out that one can derive the same condition by demanding that the OPE expansion 
should have the same form with different choices of reference square, see figure 3. It should 
have the same form, but with different OPE coefficients. 



3 Two loop polygonal Wilson loops in R 1,1 

3.1 The Octagon 

In this section we explain how to derive the octagon two loop remainder function up to an 
overall constant from the symmetries of the square together with the existence of the OPE 
limit. We start by writing down the OPE expansion of the £7(1) result (7), 



r octa 9 on = _9_ lQg (1 + e _ 2r) lQg (1 + e - 2CT) = J dpCjPip) . (9) 



The expansion (9) is a sum over twist 2k particles in the free theory where their energy 
is quantized and equal to 2k. The constants represent the amplitude for creating and 
annihilating these particles on the top and bottom parts of the octagon. They are the one 
loop form factors or structure constants. We could of course compute C^(p) by a simple 
Fourier transform 2 however we will not need to know their precise expressions. It suffices to 
notice that the k dependence is very simple, since we are expanding log(l + e~ 2r ), 3 

c£\p) = fc^cw&o (10) 

We can now move to two loops. At two loops a few things change. On the one hand the 
single particle form factors get corrected. On the other hand the single particle energies get 
an anomalous dimension of order g 2 . This second contribution gives a term linear in r of 
the form 

(r°^X aeaiiaT = -AEx 6 * / d P CM(p)e^ k (p) (11) 

k=l J 

which comes from expanding the exponential e~ T€k in powers of g 2 and using (8). According 
to (8) the anomalous dimensions jk(p) is actually independent of k at this loop order. Hence 
we have 

in T = ~Alog(l + e-*)f(a) . (12) 



2 We would get C^\p) — g 2 {— l) fe /(4fcpsinh(7rp/2)) and the integral over p in (9) should be understood 
as going slightly above the real axis. 

3 The function log(l + e~ 2r ) is nothing but the SL(2) conformal block for a dimension 1 primary. This 
is what we should expect given the discussion of the previous section. Therefore, to a certain extent we do 
not even need the U(l) result for this argument. 
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At this point, we should in principle compute the anomalous dimension j(p) and compute the 
function /(er). There is an intuitive shortcut to the correct result. The remainder function 
R, is a cyclic invariant function and therefore r <H- — r symmetric. At two loop, its term 
linear in r is given by (12). This allows us to make the following simple guess for the full r 
dependence 

Rlt%: = -y log(l + e 2T ) log(l + e~^)f{a) (13) 

Moreover, R is parity invariant and therefore r -h- o symmetric. Now we use that symmetry 
to get /(er), 

Rfi£? = -y Ml + e 2 -) log(l + e"*) log(l + e 2 ^) log(l + e" 2 ") (14) 

is agreement with [2]. The argument we presented, however, does not fix the overall multi- 
plicative constant in (14). 

To do that, we need an honest computation of ^(p). It is however easy to see that the 
normalization we picked is the correct one: Fourier-transforming f(a) from (14) we can read 
the dimension of the twist 2k fields, 

e fe (p) = 2k + 2g 2 [^(1 + ip) + ${1 - ip) - 2^(1)] , (15) 

where ip(u) = (logr)'(w). The normalization of the g 2 term in the dispersion relation is 
directly connected to the normalization of (14). Since we can compute (15) independently 
using M = 4 integrability we can fix the normalization of (14) 4 , see appendix B for details. 

It is rather simple to compute f(a) from the anomalous dimension (15): the fact that 
the anomalous dimension is diagonal in momentum space implies that in position space it 
acts as a convolution kernel. By Fourier transforming (15) we find that the logarithmic term 
of the two loop octagon is 

/^octagoBX _ r 2 r dt [2Ri(t, o) - Ri(t, a + t/2) - R^t, a - t/2)} 

V 2 loops ) linear in r T ^ gi _ \ ^ ' 

= _^ r log(l + e" 2 01og(l + e- 2CT )log(l + e 2CT ) (17) 

where R\ is the one loop answer (9). Symmetrizing we get to (14). 

These arguments allow extra terms with an OPE expansion without linear terms at large 
(positive or negative) r or a. A simple example of one such term is 

vr 2 log (1 + _L_ > ) log ( 1 



cosh 2t J \ cosh 2a 

This term respects the r — > — r, a — > —a and r -H- o symmetries of the square and has an 
OPE expansion without linear terms. Our previous argument would be blind to it, we would 



4 A shortcut is to recall that the large p behavior of the one loop anomalous dimension of excitations 
around the GKP solution should diverge as Tlog(p) where T is the cusp anomalous dimension which is 
known [7, 8]. 
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Figure 4: The Decagon null Wilson loop in blue embedded in the Penrose diagram of R 1 ' 1 . 
The red dashed line is the reference square. 



have to resort to a physical argument to discard such term. In this particular case we could 
discard it since it has extra singularities when r = ±7r/4, ±37r/4 etc, which should not be 
there. It would be interesting to understand in more generality to which extent is M = 4 
OPE constructible. 

In fact, the OPE is very similar to looking at particular physical cut contributions in 
amplitudes. It would be interesting to find a precise relation. 



3.2 The Decagon 

In this section we use the same technique to derive the Decagon two loops remainder function. 
The generalization to any even n > 10 is given in Appendix A. We start by gauge fixing the 
decagon cusps to be at 

x+ = {-oo,- X + --^,0,1} (19) 
- { 1 + X2 1 n 11 

x i ={-007 — , — -,o,i}. 

Xi Xi 
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The cusps are located as in (5). A corresponding basis of four conformal cross ratios is 



Xi 



X 5,l X 4,2 
X 4,5 X 1,2 



Xi 



X 4,1 X 3,2 
X 3,4 X 1,2 



Xl 



X 1,3 X 5,A 
X 5,1 X 3,4 



% A 1 OC 



x 2 



4,l lX '3,2 



X Q A X 



(20) 



3,4**'l,2 



Any OPE limit corresponds to a choice of two non adjacent edges in the x + or x~ directions. 
For the Decagon there are ten different such choices, i.e. ten different OPE limits. These 
limits are all related by the symmetries of the Decagon which reshuffle the cross ratios. These 
symmetries are parity, which reflects the sign of time, and cyclicity, which cyclicly permutes 
the labeling of the x + and x~ points. Combining these two we get a Zio transformation that 
acts on the cusps points as 

(21) 



Xi — > x 



i+1 ' 



Under this transformation we have 5 
1 + X 2 " 



Xi -> 



Xi 



Xi 



1 5 



xj — > x. 



Xi -> X 2 



X 2 -> 



i + xi 

Xi 



(23) 



We will now construct two loops Decagon OPE expansion in one channel. Demanding 
the answer to be Zio symmetric will guaranty the correct OPE limit in all other channels. 

The U(l) part. We choose the two null lines for the OPE expansion to be in the x 
directions and located at x + = and x + = 00. As before, we choose the reference square 
cusps on these lines to be located at x~ = and x~ = 00 (see figure 4). For that choice, 
the OPE expansion parameter is Xi = e ~ 21 '■> whereas xt = e~ 2cr . The U(l) result in that 
channel is read from the one loop amplitude as explained in section 2. It is given by 



rum 



log 



1 + 



Xi 



+ ylog 



1 + 



1 + X 2 ~ 
Xi 

1 + X2 J 
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log [1 + Xi]-~2 lo S i 1 + Xi] lo S 



1 + 



Xi 



Xi 



(24) 



log 



1 + 



Xi 



+ [terms independent of r] 



1 + 

To see that, note that the Decagon can be decomposed into three Octagons as in the bottom 
of figure 4. Since the U(l) part is linear, it decomposes into a sum over these three Octagons 
corresponding to the three terms in (24). The same applies to any n > 8 (see Appendix A). 

The OPE procedure. The expression (24) is written as a sum over three terms, each 
of the form of the Octagon U(l) part with shifted r and a. That is, each term has the form 

dp e ~ ip{a+b/2) 

k " 



log (1 + e~ 2T+a ) log (1 + e- 2 ° +b ) = J2 {-e- 2T+a ) 



k=l 



2(p-z0)sinh(7rp/2-z0) 



(25) 



5 The best way to find how that transformation acts on the conformal cross ratios (20) is to map them to 
the Y functions of [9] . The reason is that (21) has a simple action on these cross ratios given by Fj a ' — > Y"J a+1 ' . 
For the Decagon Yq = Y3 = and yj a+10 l — Y~J a '. The translation between the Y functions and the cross 
ratios (20) can be read from [9] and is given by 



+ _ r [i] _ Y [t 



v+ — v — vl 5 l 



y~ - Y, - Y [5] 

Al 1 1 1 2 > 



- yj 11 = y! 6] 



x 2 =H' = H" ■ ( 22 ) 

We can then use the Y system equations Y"J a+1 'yJ a 1 ^ = (1 + Y"j"\)(l + vM) to map all the Y functions 
generated by the Z 10 action to the ones appearing in (22). Using these equations one finds that (21) acts on 
the cross ratios (20) as (23). This procedure generalizes immediately for larger number of edges. 
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As for the Octagon, we identify the terms e~ 2kr in the sum with the exchange of the twist 
2k one particle state. At two loops, that term will give rise to a term linear in r obtained 
by dressing the previous expression by 



-T7i(p) 



-2Tg 2 [4){l + ip/2)+4){l-ij)/2)-24){l)} . (26) 
log (1 + e~ 2T+a ) log (1 + e~ 2(7+b ) ^ 2/rlog (l + e~ 2T+a ) log (l + e~ 2a+b ) log (l + e 2a - b )(27) 



We conclude that 



By appling (27) to the three terms in the U(l) part (24), we find the part of the two loop 
remainder function which is linear in r: 



d(2) 
linear in r 



log [xi ] log 



-y log [Xi] log 



+ y log [xr] log [! + Xi] log 



1 + 



1 + 



Xi 



I + X2 

1+X2 



log [l + Xi"] log 

xi 



1 + 



1 



(28) 



log 



1 + Xj 



xr 



X2 



log 



log 



1 + 



1 + X2 + 

xi 



Xi 



Xi 



Symmetrization. This expression is not invariant under the Z10 symmetry and therefore 
cannot be the full result. We should now proceed like in the octagon case. I.e. we should first 
replace log(xr) by some function which behaves as log(Xf) when \i ~* 0- Next we should 
symmetrize the result over the Z 10 action. Now, in general, this symmetrization will spoil the 
OPE expansion (28). The idea of the bootstrap program is to look for a replacement which 
preserves (28). Remarkably, the simple replacement \og{x\) — ► — l°g(l + VXi ) °- oes * ne J 00 ' 
Of course, by construction, the remainder function obtained in this way is Z10 symmetric. 
We therefore arrive at the following expression for the two loop Decagon remainder function 6 



R 



2 loop 



-ylog 



9 1 

-ylog 



-ylog 



9 1 

-ylog 



9 1 

-ylog 



X 2 



Lxi 
1 

LX2 
1 

xr 
1 

xr 
1 

X2 



+ 1 
+ 1 

+ 1 

+ 1 

+ 1 



+ 1 


log 


f xr +1" 


log 






.xr + 1 . 





(X^ + 1)(X 2 + + 1) 



xixi 



log 



( X + + 1)( X2 + + 1) - 
Xi" + Xt + 1 . 



log 
log 
log 
log 



(xr + 1)(X 2 +1) 

xr + xr + 1 

(xr + i)(x 2 +1) 
xr + xr + 1 



log 
log 



xr + i 



1 



X 2 

xi + xi + i 



xi 



X 2 



+ 1 



log [xi + l] log 
log [xt + !] log 



xr 
xr 



log 
log 



xr 



Lxr + i 

xi 



+ 1 



(29) 



xr + i 



+ 1 



x 2 +1 

xr + x 2 + 1 
xr + 1 



which is compatible with the OPE expansion in all possible channels. Equivalently, 



R 



2 loop 



9 H 



log 



1 " 

1 + — 


log 


[l+ Xr 1 


Xi . 




. I + X2". 



log [l + xi] log 



1 + 



+ cyclic permutations 



J This result was independently derived before [4] appeared. 
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where cyclic permutation stands for the Z 5 C Z 10 cyclic permutation of the x + and x~ cusps. 
Up to a — 7r 4 /12 constant, that is indeed the Decagon two loops remainder function guessed 
(and confirmed by numerics) in [4]. In Appendix A we generalize this result for n > 10. 

3.3 General gauge theories 

In the above results we have only used the symmetries of the problem. Thus, we expect 
such results to be valid in any conformal gauge theory with a weak coupling limit. It is also 
easy to understand from diagramatic point of view why the results are valid for any planar 
conformal gauge theory at two loops. Consider all the two loop diagrams. Many of the 
diagrams involve only gauge bosons and are the same in all theories. The one diagram that 
is different is the bubble correction to the propagator. Such a diagram contains a piece with 
gauge bosons and a piece that involves a matter loop. The theory is a gauge theory coupled 
to matter in such a way that it leads to a CFT with a tunable coupling, with zero beta 
function (by assumption). The matter theory has an global symmetry that we are gauging. 
Let us call j the corresponding currents. We have a current two point function (jj) ~ k/r 6 . 
The condition that the /3 function vanishes is the condition that the constant k has the value 
it has for the matter in M = 4 SYM. Thus the diagrams involving the matter bubble have 
the same value that they would have in M = 4 SYM. This argument also appears to work 
for non-planar gauge theories (at two loops). In this case we have different diagrams, but 
only involving gauge fields. 

It would be interesting to investigate a similar argument for the case of three dimensional 
conformal field theories, such as ABJM [10]. 

4 Remarks on the three loop octagon 

We can attempt to extend this method to the computation of the octagon Wilson loop at 
three loops, still in R 1,1 kinematics. One difficulty is that at two loops we might be able to 
create propagating states consisting of two particles. Let us ignore this issue and let us first 
focus on one part of the answer which only involves the propagation of one particle states. 

At two loops we obtained the part of the two loop answer that goes like tF 2 (t, a), where 
F 2 has an expansion in powers of e _2r . This was computed by convolving the anomalous 
dimension kernel with the one loop result, see (16). At three loops we can easily obtain 
the piece that goes like t 2 F^{t, a). We simply need to convolve twice with the anomalous 
dimension kernel appearing in (16). We obtain 



Even though it is not manifest, we have that /3(a) = f^{—a). If we now symmetrize this 



r 2 F 3 (r,a) 



^7 * 7 * Ri = y (2s 2 )V Iog(l + e~^)h{a) 
-2Li 3 (-e 2CT ) + 2<xLi 2 (-e 2CT ) - ^ log 3 (l + e 2CT ) + 
+4alog 2 (l + e 2CT ) - (2a 2 + ^) log(l + e 2CT ) 




(30) 
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under r -H- cr we obtain 

^ oc // 3 (r)/ 3 (a) (31) 

As we take r — > oc we get the right r 2 terms, in agreement with (30). One could wonder 
if (31) is the full answer or not. If one assumes that the three loop answer factorizes into 
a function of r and a function of a, then (31) is the only consistent answer. In fact, if we 
assume that the answer factorizes at all loop orders, then one can determine it using the 
same method. However, we do not know any reason why it should factorize. And the strong 
coupling answer [11] appears inconsistent with this factorization hypothesis. This does not 
say at what loop order factorization would stop. There is another potential problem with 
(31). We can consider the term linear in r. The guess (31) gives a very specific term linear 
in r. On the other hand the OPE analysis gives two sources for terms linear in r. We 
have one loop corrections to the OPE coefficients and two loop corrections of the anomalous 
dimensions. The latter on its own gives a term which has a rather different structure, being 
of transcendentality two in r and four in a. However, we have not computed in detail the 
corrections to the OPE coefficients. This correction could involve the creation of two particle 
states and we have not fully analyzed it. We leave this to the future. The conclusion is that 
probably (31) is just a part of the three loop answer, which correctly captures the full r 2 
behavior, but by some miraculous cancelations it could be the full answer. 



5 Conclusions and open problems 

In this paper we have provided a quick derivation of the two loop result for polygonal Wilson 
loops in R 1 ' 1 kinematics. It is based on an important feature of the OPE expansion of the 
two loop results: the terms which have a logarithmic discontinuity around the factorization 
locus can be predicted from the OPE expansion of the one loop result. The idea is to start 
from the known result for the one loop Wilson loops. This is a simple function which is 
just the result that we would obtain in a U(l) theory. We then apply the Operator Product 
Expansion for Wilson loops derived in [1]. One selects a particular expansion channel and 
the structure of the OPE, together with the knowledge of the anomalous dimensions, allows 
us to derive certain logarithmic terms in the two loop result. They are the discontinuity 
of the two loop answer around the factorization locus. One then guesses a simple way to 
complete these terms into a full expression. When we select an expansion channel we break 
the cyclic symmetry. Imposing the cyclic symmetry gives us a natural way to complete the 
answer into the full expression. In this way we rederived the two loop octagon expression 
(14), derived originally in [2]. We have extended these results for higher number of gluons. 
These expressions were guessed previously in [4]. Their approach was different. They just 
assumed a form for the functions that could appear in the answer and then imposed the 
correct colinear limits, which would be the zeroth order term in the OPE. In fact, if one 
knows the functions that appear in the answer, then one can find the right combination by 
looking at simple limits. On the other hand, the OPE, can be good for determining, or at 
least constraining the functions that appear. In this case with R 1,1 kinematics the OPE gives 
us the functions that appear in the answer. 

It is natural to try to extend this two loop analysis to the full R ' 3 kinematics. In principle, 
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the same method works. In practice, it is more complicated because there is a large number 
of SL(2) primaries that appears. In the case of R 1 ' kinematics, we had simply one tower 
of primaries which are simply plus and minus derivatives of In the R 1 ' 3 case we have 

primaries of the form D l z F +z and D l s F +z -, where z and z are the transverse coordinates. The 
I index is related to powers of one of the extra angular coordinates that we have when we 
do the OPE in the R l,a context. However, we also have primaries of the form D l z F_ z and 
DlF_g. When we perform the expansion of the one loop answer, we need to separate the 
contribution from these two sets of primaries. We leave this to the future. 

We can also consider higher loop contributions. The leading logarithmic terms are easy 
to obtain, one simply repeats the convolution with the one loop anomalous dimension kernel 
once again. For R 1 ' 1 kinematics we wrote an expression which reproduced all the r 2 (and 
a 2 ) terms in the OPE at three loops. It remains to be seen if there are other terms that we 
should add to that expression in order to reproduce the full answer. In principle, we can also 
get subleading logarithmic terms. These arise from higher loop corrections to the anomalous 
dimensions, which we can compute, and also from corrections to the OPE coefficients. These 
are harder to compute since they might involve multiparticle states, etc. In short, it would 
be great to gain analytic control over the breaking of the SL(2) symmetry discussed in the 
text. This seems to be the main missing link to be able to extend our techniques to any 
loop order. Of course, if one restricts to the lowest terms in the OPE, which only receive 
contributions from single particle exchanges, then it is easy to go to higher loops once one 
knows the lower loop answers, as explained in [1]. 

Recently the full TV = 4 all loop integrand was proposed in a remarkable Yangian invariant 
form [12]. It would be interesting to see if the OPE limit can be applied directly at this level. 
Also, it was recently understood how to generalize the bosonic null Wilson loops into super 
loops which are dual to amplitudes with arbitrary polatizations [13, 14]. It would be very 
interesting to apply the method we described to compute non MHV scattering amplitudes. 
In these cases, fermionc excitations of the flux tube will also be excited. Another connection 
which would be interesting to work out concerns the relation between the OPE expansion 
and the high energy Reggee limit of scattering amplitudes. They seem to be closely related. 
The recent papers [15] seem to be a promising starting point for establishing a more precise 
connection. 
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Figure 5: Null Wilson loop with n = 2m edges in blue embedded in the Penrose diagram of 
R 1 ' 1 . The red dashed line is the reference square. 

A Polygon Wilson loops with n > 10 sides 

In this Appendix we use the consistency of the OPE expansion to obtain the two loop 
remainder function for any null polygon in R 1 ' 1 kinematics. As for the octagon and decagon, 
we start from the U(l) contribution and use the knowledge of the one loop energy (15) to 
obtain the terms linear in r in all channels. We will then show that the guess of [4] is the 
simplest solution that is consistent with all the OPE expansion in all channels. 

Consider a n = 2m null polygon in R 1 ' 1 . An OPE channel in the x + direction is defined 
by splitting the cusps points xf into two groups of ordered points, top and bottom. We denote 
these by (xf, x^ +1 , . . . , xf ) and . . . , x+) respectively, see figure 5. The locations of top 

and bottom cusps in the x~ direction are (x~ , . . . , xj) and (x~ , . . . , We chose the 
reference square cusps to be at (x^,xj), and (x^xj). 

To obtain the U(l) contribution in this channel we decompose the n-gon contour as a 
sum of octagons as we did at the bottom of figure 4. The cusp of the octagons are denoted 
by A, . . . , H, see figure 6 for an example. In this decomposition all octagons share the two 
cusps of the reference square A = (x+ , x~ ) and E = x~). Each octagon also shares one 
more cusp from the top group 

C G {(x 1 , X t ), , X l ), (X2, X 2 ), • • • , (^ s _i) ^5-2)5 x s-l)} (32) 
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Figure 6: Example of a Null Wilson loop (thick dark blue line) together with a reference 
square (dashed red line) and one of the octagons (thin light blue line). The points A and E 
are cusps of the reference square. Points C and G are points of the top and bottom part of 
the polygon respectively. As depicted in this figure, the points B, D, F, H are automatically 
fixed given these four points. 



and one more cusp from the bottom group 

G ^ {( X 8+V X s+l)i ( X t+2> X s+l)' ■ • • ' ( X m-1? X m-l)} (33) 

So far we described four of the eight cusps of each octagon. It is easy to see that the four 
remaining cusps are uniquely fixed once these points are chosen, see figure 6 for illustration. 
Although we will not explicitly use them in what follows, their positions are 

B = (C + ,A~), D = (E+,C-), F=(G + ,E~), H=(A + ,G-). (34) 

We have therefore a unique decomposition given by a sum over the points C and G (all 
other points are fixed in terms of these). The sign of each octagon is determined from its 
orientation and is given by (— l)( c+ ~ c ) where c + stands for the index appearing in 

C + , i.e. C + = etc. To summarize, the U(l) part in that channel is given by 

<?) = ~^- i y + ^ +g+ ~ 9 ' l ^( 1 + Y cAEc)^g(l + Y^ EAC ) . (35) 

C,G 

where Y ABCD = [^±Ip±j[^±Ig±j ■ Note that Yq AEC oc e~ 2r and Yqeac ^ e ~ 2a - We can now 
use (27) to derive the term in the two loops remainder function that is linear in r 

in r = 9*t J2hl) c+ - c ~ +9+ - 9 ~ log (1 + Y+ AEC ) log (1 + Yo SAO ) log (1 + l/Ya EAC ) . 

C,G 

(36) 
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The full two loops remainder function is a function whose term linear in r in any channel is 
given by (36). To guess this function it is useful to notice that the sum over C and G can 
be partially performed. First we notice that 

1 + Y GAEC = 1 + / + _ + \/ + _ + \ = ( + _ T+ U T + _ „+\ = 1+ ~+ \*') 
\ U/ g+ x s )\ Jj c+ ^m) \ g+ s )\ x c+ ^m) X g+,s Jj c+,m 

Next notice that, from (32), we have Ylci- • • ) = J2 c +=c- c -+i (■ ■ •) anc ^ similELX" for G. 

We should be slightly more careful with the boundary terms, e.g. when c~ = s — 1 we only 
have one allowed value for c + , namely c + = s — 1, see (32). In other words, the sums over 
c + and g + are rather trivial as they contain either two or one terms. E.g., 

\^ \^ (-iy + -c-+ S + - a ~ log V.^,™ = log V + i,c- + iV,c- | (38) 



g + ,s c + ,m I \ 9 +l) c 9 +l,c 



which is valid away from the boundaries, i.e. for c < m — l and g < s — l. At the 
boundaries, if g~ = m — l we have 

£ (-1)^-) log = g + " +1 ' C ' +1 + g '- C ' (39) 

z — ' \ X + X + / \ X _ , -, _x _ t1 _ / 

c+=c- \ 9 + . s c+,m/ \ g +l,c g +l,c / 

with a similar expression for c~ = s — 1. Results (38) and (39) are quite interesting in the 
sense that they involve the logs of nearest-neighborhood cross-ratios, the indices s and m 
dropped out completely. These building blocks posses a very useful property: 



log 



X a+l,b+l X a,< 
X a+l,b X a,b+l 



gives rise to a linear r behavior if and only if a = s and b = m or viceversa. We can use this 
to write a cyclically symmetric sum which reproduces the r- linear part in the large r limit: 

ST^ j Q [ X t+l,c+l x a,c \ j q / x b+l,d+l X b,d \ j q | 1 + ^M^V J J Q [ ]_ _|_ X a ^ X d > c 

a,6,c,d in cyclic order \ a+1 > c a ' c+1 / \ X b+l,d X M+l/ \ X a,b X c,dJ \ X a,d X b, C/ 

Although it is not manifest, this expression is symmetric in exchanging x + and x~ . Indeed, 
it is simply the result of a telescopic sum over x~ of the expression guessed in [4], thus 
showing that it satisfies the r-linear part of the OPE requirement in all possible channels 
and is therefore the simplest and most natural solution. 



B Excitations from Integrability 

In this appendix we compute the anomalous dimensions of the flux tube excitations using 
integrability. More precisely we will use the mapping of the flux tube to the GKP state 



ZD+...D+Z + ... (40) 
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and study the energy of fluctuations around that state. Here D + is a light-cone direction 
and Z is one of the complex scalars of the theory which need to be included so that the 
derivatives have something to act on. We can think of the scalars as being two fast particles 
sourcing the flux tube represented by the light-cone derivatives. The spin S is to be taken to 
infinity as explained shortly. This formula is of course schematic, the state is a complicated 
superposition as indicated by the dots and an overall trace is omitted. 

We are interested in excitations 

ZD + ... X ...D + Z + ... (41) 

moving in the background of derivatives. The effective length perceived by these excitations 
is 2 log S. The energy of this state is given by 

E = 2 lcusp \ogS + C + Yl e x(Px) ( 42 ) 

excitations \ 

where p x is the momentum of the excitation and e x (p) is the dispersion relation of the 
excitation \- The first term represents the vacuum energy contribution and C is an irrelevant 
constant, also associated to the vacuum. There are several different excitations we could 
consider. In this appendix we will consider scalars, fermions and, most relevant for this 
paper, the case when 

X = D k _F + - . (43) 

This is the most relevant case since excitations of the form (43) are the ones that we expect 
to generate in the OPE Wilson loop when considering loops in R 1,1 as in this paper. Of 
course the detailed structure of the Bethe state with excitations (43) is quite complicated 
since it will mix with fermions and scalars. For us what is important is that is has the correct 
quantum numbers to be created by the polygon Wilson loops. For the bulk of the text, the 
main result is that these excitations have a one loop energy given by 

tD*F + - (p) = 2 + 2k + 2g 2 (V (l + *!) + V (l - i|) - 2^(1)) + 0{g A ) . (44) 

Excitations of the GKP string were considered in much greater detail by Benjamin Basso 
in [16] who provided all loop expressions for most excitations. In this appendix we will 
consider mostly one loop and will mainly focus on the excitations needed for the main text. 
We will also use the method of Baxter polynomials which is very suitable for perturbative 
weak coupling computations. 

B.l Scalars Z at one loop 

To construct the state with a single scalar x — Z moving in the see of derivatives as in (41) 
we use Beisert-Staudacher (BS) Bethe equations [17]. For completeness these equations are 
written down in section B.5. The precise notation is introduced in that section. 

These equations describe excitations around the Z vacuum. From this point of view the 
excitations are the light-cone derivatives while the vacuum is made out of scalars. In other 
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words, (41) is an extremely excited state. Luckily there is a way of efficiently exchanging 
particles and vacuum constituents in Integrable models via what is called a particle-hole 
transformation. This can be nicely implemented using the Baxter polynomial approach as 
proposed in the AdS/CFT context in [18]. We will follow this approach in what follows. 

To build (41) we consider BS equations with S momentum carrying excitations (the 
derivatives) and L = 3 (the scalars). At one loop the relevant Bethe equations are equivalent 
to the Baxter TQ relation 

T(u)Qi{u) = $(u + i/2)Q 4 (u + i) + $(u - i/2)Q 4 (u - i) (45) 
where $(w) = u L with L = 3, Qa{u) = Ylk=i( u — u^) and 

Kg 

r(«) = 2n(t*-^). (46) 

3=1 

The zeros of T(u) are called holes. The index 4 comes from a conventional labeling of the 
Bethe roots in the BS equations, see section B.5. Since T(u) is a polynomial, the left hand 
side of (45) vanishes for u = which coincides precisely with Beisert-Staudacher equations 
when only Z's and Z) + 's are present in the single trace operator. Note that for (45) to make 
sense we must have 

K e = 3. (47) 

The idea is to trade the dynamics of the S particles Uj by the dynamics of the 3 holes 9j. 
Actually, as explained in section B.4, two of these holes are located at 

B*±S= (48) 

and do not have therefore any interesting dynamics. They are the so called large holes. 
This result is valid for any coupling, with any finite number of excitations with arbitrary 
polarizations. The remaining hole is the interesting one, the small hole. It corresponds to 
the excitation x = Z moving in the vacuum. The goal is then to find the momentum and 
energy of this holes, pz(0) and €z(Q) to derive e^(p). 

At values of u with positive imaginary part the first term in the right hand side of (45) is 
way bigger than the second one which can therefore be dropped'. The resulting functional 
equation is then easy to solve [18], 

3 r 



QAu) = c2 u l % TT — j \- , for u in the upper half plane , (49) 

AJ - p ( u+i/2 \ 

3=1 L { * ) 

where c is a constant which is irrelevant for our argument. For negative imaginary values 
of u we would drop instead the first term in the right hand side of (45) and we would get 
the complex conjugate function. We can now compute the energy and momentum of the 
excitation. 

7 This is so for |u| <C S which is more than enough for our purposes 
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The energy of the state is given by (74). When taking the logarithm of Q4, the product 
in (49) becomes a sum of three terms. Two of them are associated to the large holes (48) 
and gives us the vacuum energy contribution, containing a term proportional to log S plus a 
constant. The remaining term gives us the energy of the excitation. More precisely, we find 
(42) with lcusp = Ag 2 + 0(g A ) and [18] 

e z (0) = 2< 7 2 (V<l/2 + i0) + ^{1/2 -id) -2^(1)) + 0{g 4 ) , (50) 

The momentum is read off from the effective Bethe equation for the hole excitations. These 
are nothing but 



= — 1 , at u=6L- (51) 

where f ± = f(u ± i/2), f ++ = f(u + i) etc. This is of course the definition of the holes. 
When plugging (49) into this equation, the left hand side will contain a factor of the form 
e -ip z (9)x2iogS g enera ted by the large holes in (49). From this term we read off the momentum 
Pz(0) (recall that the effective length of the operator is 2 log 5). In other words, 

Pz{0) = hm — log — — . 

5^00 2 log S Qa\v—i) 

We find 

Pz (6) = 29 + 0{g 2 ). (52) 
Combining this result with (53) we arrive at the known result 

e z (p) = 2g 2 &(l/2 + ip/2)+ ^(1/2 - ip/2) - 2^(1)) + 0(g 4 ) . (53) 



B.2 Other excitations 

To consider other excitations in Beisert-Staudacher equations associated with the Dynkin 
diagram of figure 9, we also need to consider the auxiliary roots Uj for a / 4. The Baxter 
equation (45) for the momentum carrying roots is now replaced by 

TQ 4 = $ + Q++Q 3 -Q 5 - + Q-QrQtQt (54) 

where T{u) is given by (46) with Kg = L + K3 + K5. The excitations and are 
fermionic excitations. Note that when u is in the upper half plane the first term in (54) 
dominates exponentially over the second one. However, when u ~ + i/2 with a = 3, 5 
this term is suppressed by Q~ . This means that, with exponential precision, there is a 
complex zero of T{u) for each fermionic root located at +i/2. Similarly, there is another 
hole at — i/2. 

There is an important exception: the case when two roots of type and u^p are in 
the same position v. In this case we also have the possibility of having one hole at v + i/2 
and one hole at v — i/2 instead of two holes at v + i/2 and two holes v — i/2, as one can 
conclude after a simple inspection of (54). Hence, from now on, we shall use Vj to denote the 
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roots u~p = in this particular case, i.e. the which are at the same position and which 

are accompanied by only two holes at Vj ± i/2. At the same time we shall use Uj and u^p 
to denote the remaining roots of type 3 and 5. Then, 

T(u) = 2 U - ^ ) J] ({u - zf + -\ H(u-z) (55) 

J z6{ Vj }U{«f }U{uf } 2 ={^} 

where (9j are real small holes like the one studied in the previous section which are not bound 
to any fermionic roots. Now, as before, we can solve the functional equation (54) by going 
to the upper half plane and using simple functional identities of Gamma functions. We find 

1 q + 

Q±{u) = c exp (— 2ui log S) — q$ <lt it > f° r u m ^ ne u PP er half plane (56) 

where = Y{u/i) L and 

/ u _J-)\ 

9a = TT r — — • (57) 




Note that Q\ vanishes for u = v + i/2 which means that there is a complex root at this 
position. Similarly from the expression for Qa{u) in the lower half plane we conclude that 
there is another root w 4 at v — i/2. Figure 7 summarizes the three possibilities we found so 
far: single real hole, fermionic excitation dressed by two complex holes and double fermionic 
composite dressed by two complex holes and two complex momentum carrying roots. In this 
picture we only represent the holes and the particles u^ a ' with a = 3, 4, 5. All other auxiliary 
roots don't couple directly to the momentum carrying roots and hence carry no energy 
nor momentum at this loop order. When considering two loops and higher we would need 
to consider the energy of excitations u^ 1 ' and as well. 

We can now use again (74) to read of the energy in the presence of all possible excitations. 
We find 

E - E vacuum = 2j e i>( z ) + e ^( 2; )' ( 58 ) 



where (53) and the energy of and excitations are given by 

e^fl) = 2/(^(1 + ^)+^(l-^)- 2^(1)) + 0{g A ). (59) 

Note that the excitations Vj carry no energy at all!, they are not present in (58). We can 
now compute the momentum of the several excitations. The momentum is read off from the 
term e ip ^ 2logS arising in the corresponding Bethe equation in (75). For composite states 
we should multiply the corresponding equations as usual. Plugging the corresponding Q 4 
given by (56) in (75) we read the momentum of these fluctuations as before. For and 
excitations we find again 

p46) = 26 + 0(g 2 ) (60) 
while the v excitations have carry no momentum, 

p v (6) = + 0{g 2 ). (61) 
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Figure 7: Fundamental building blocks using the holes and the fermionic roots and . White 
balls represent holes; Blue balls correspond to fermionic roots ; Red balls indicate fermionic roots 
M ( 5 )j Green balls denote complex momentum carrying roots . The energy of the first building 
block (the single hole) is given by (53) while the energy of the fermionic excitations in the middle is 
given by (62). The last composite, to the right, is isotopic to all loop orders [16] and hence carries 
no energy or momentum. 

In sum we found 

e^(p) = 2g 2 [^(1 + ip/2) + ^(1 - ip/2) - 2^(1)] . (62) 

Using these results we can derive the dispersion relation of several different excitations. First, 
consider the excitation associated with a single root #. An excitation of this type carries 
the quantum numbers of C4&1, see figure 9. This bilinear converts 

D + Z .D 1:L $43 <-> fclatcJcJlO) (63) 
into a twist one fermionic excitation 

*i 3 «S4lo) 

Hence we derived that 

e^twist onc(p) = 2g 2 [i/>(l + ip/2) + V(l - ip/2) - 2^(1)] . (64) 

Suppose we added a bosonic root to this excitation. This would amount to acting with 
the oscillators Csc\ which would convert the R-charge index 3 of the fermion into a 4. This 
would be another fermionic twist one excitation. Since excitations do not carry energy 
or momentum this excitation would have the same anomalous dimension which is of course 
what we expect. 

Lets now consider a real hole 9 together with a composite of type v described above. As 
we saw the v composite carries no momentum or energy at all. There is a simple reason [16] 8 
for this which becomes clear when we identify the excitation using the oscillator picture. By 
adding one root and one root we convert the scalar excitation Z = $43 into another 

8 We thank Benjamin Basso for illuminating discussions on this point. See [16] for more details. 
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Figure 8: Pattern of Bethe roots corresponding to the excitation D_J-aa. This pattern of roots is 
what we would get if we would add small twists to the Bethe equations to regulate them. Otherwise 
we would find instead a descendent of D^^Faa. In other words, some of the roots in this figure would 



go to infinity. The excitation D k _T\2 is trivially related to this one by a wing exchange Uj <->■ u 
It has therefore the exact same energy and momentum. The state (43) is nothing but the linear 
combination D^T\2 + D_J-a 2 an d therefore it has the same energy as either of these states. Now, 
at one loop, only the first constituent in this figure - the single root it' 3 ) - carries energy and 
momentum! All other excitations are isotopic. Hence (44) follows. 



,(8-a) 



scalar excitation Y = $32. The v composite is an isotopic degree of freedom responsible for 
implementing the 5*0(6) symmetry of the GKP vacuum! This is why it has zero energy and 
momentum as we checked at one loop order. Of course, from this symmetry argument the 
v composite ought to have zero energy and momentum at any loop order and it does indeed 
[16]. ' 

Finally we could consider more complicateds excitation such as (43). It is described by 
several roots of all possible kinds, see figure 8. However, we see that only a single fermionic 

root carries energy at one loop order. Hence (44) follows. At higher loops the roots u^p and 

(7) 

are no longer isotopic in this vacuum. Thus the state (43) no longer has a reasonable 
description in terms of a single particle state; instead it decays into several excitations (e.g. 
fermions). 

Now, this result is part of a much more general observation which is that at one loop the 
roots and vP' are isotopic and hence there is an isotopic symmetry enhancement which 
will be broken at higher loops when these particles become momentum and energy carrying. 
This symmetry contains the SL(2, R) symmetry mentioned in the main text. It implies that 
the one loop spectrum is organized according to SL(2, R) primaries and descendents. For 
example, all excitations D^_ip (or D^X) have the same anomalous dimension because the 
pattern of Bethe roots associated to this excitation is just the Bethe roots of ip (or X) plus 
a bunch of isotopic roots as in figure 8. At higher loops the roots vS 1 ' and start carrying 
energy and momentum and all these degeneracies are lifted. 
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One loop Summary 



We have found that the one loop dispersion relation for excitations around the flux tube is 
always of the form 

e s (p) = 2g 2 [i>(s/2 + ip/2) + ^(s/2 - ip/2) - 2^(1)] , 
where s = S + A is the conformal spin of the excitation. For example 

s excitations 
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1p2a 1p 2 a 
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B.3 Higher loops 

In this section we review how to set up the computation of the anomalous dimensions at 
higher loop orders. We will illustrate how to perform this computation at 2 loop order but 
the generalization to higher loops can be done along the same lines. For all loop results see 
[16]. Up to four loops the Baxter equation reads 



® + Qi + QzQi 



{B)? )+ yBtBtBtB$ 



+ c.c. 



(65) 



See section B.5 for the definition of the several Baxter polynomials. Up to wrapping order 
g 2L , the right hand side of this equation is a polynomial in u which allows us to define a 
polynomial T(u) thus simplifying the computations. This justifies the splitting of terms 
between the displayed terms and those in "c.c". Now we proceed as before. For a very 
similar computation see [19]. For positive imaginary part of u we can drop the "c.c." terms 
in (65) and find 



Qi{u) 
Q 



(0), 



l + g 2 i 



u 



Ld u + 2d u \ogQf 



u=i/2 



+ ^2d u log Q a 



a odd 



u=0 



Ml - in) + 0(g A ) , (66) 



where Q^\u) is given by (56). As before, the energy of the several excitations can then be 
extracted from (74) while their momentum is identified by looking at the term proportional 
to 2 log S in the (logarithm of) the corresponding Bethe equation. For holes and fermionic 
excitations of type and the formulae of the previous section are corrected to 



where 



e z (u) = + lcusp ^ i h \l/2 + iu) - 2g 4 4 + \l/2 + tu), 
ty(u) = +lcus P ^ { o + \l + iu) - 2g 4 ^ 2 +) {l + iu), 

tPi + \b + iu) = \ [ip a (b + iu) + ^J a (b - in) - 2(a + l)Va(l)] 



(67) 



(68) 
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Roots of type 1 and 7 only acquire energy at 3 loops, at 2 loops we still have ei/jiu) = 0. 
For the momentum we find 



Pz(u) 

Pl/7(u) 




(69) 



B.4 The two large holes 



In this appendix we explain why there are two (and only two) large holes located at ±5/ y2. 
The argument relies on two simple observations. 



1. Consider 



^ / \^ / v^7 \ = 2cosp(w) 



(70) 



at large values of u ~ S. For example, at one loop we have, see (54), 



p{u) 




(3) 



1 KS 

-Y 

2 ^ 

.7 = 1 



1 



(5) 



(71) 



The precise form of this expression is not very relevant (e.g. at higher loops this 
expression is modified but its form as a sum of poles still holds and the argument that 
follows goes through untouched). 

For u ~ S the roots are well described by a density and the function p(u) develops 
two large cuts (one for positive u and another for negative u). The real part of p(u) 
on these curs coincides with the (logarithm of the) Bethe equations in the scaling limit 
and is therefore equal to the mode number of the corresponding roots. For the GKP 
ground state, it equals n (— ir) for positive (negative) u. For even larger w's, i.e. after 
the cut, the function p(u) decays to zero. Between p(u) = tt and p(u) — > there is 
only one value where cosp(w) vanishes which is of course p(u) = tt/2. Thus there is 
only one large hole for positive u. Similarly there is only one large hole for negative u. 

2. At large u we have 



This follows trivially from expanding (54) at large u. Now, the existence of only two 
large holes proved above implies that t 2 k/u 2k <C ^/w 2 for k > 1. Thus, for large u, 
we have T(u) oc 2m 2 — S 2 so that the two large holes, which are the two large zeros of 
T(u), are indeed located at ±S/\/2. Q 

9 It is also fun to note that at the end of the cut we have p(u) = tt which means 2 + H ~ — 2 yielding 
u ~ S/2 which is indeed the known value for the end of the large cuts. 




(72) 



where 



h = S 2 + O(S) . 



(73) 
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Figure 9: Beisert-Staudacher Dynkin diagram and corresponding oscillators. 

The Af = 4 Asymptotic Bethe Ansatz 



In M = 4 SYM single trace words are made out of letters which can be written using bilinears 
of oscillators such as 

®ab ^ 4 C t|0) 



o 4&j,|0) , ^ 6t 6 t| ) j jr. . ^ a t a t | 0) ; ^ ^ c t 6 t | ) ,etc. 



The oscillators can have SU (2) x SU(2) Lorentz indices (a a and 4) or (4) i?-charge index 
(4). The oscillators a and 6 are bosonic while c is fermionic. Difference single traces are 
made out of different numbers of bilinears, denoted by K a , see figure 9 and [17, 20, 21] for 
more details. The anomalous dimensions of states for fixed K^s are then given by solving 
K i + . . . Kj algebraic equations for the so called Bethe roots {Uj } where a = 1,...,7 
and j = 1, . . . , K a . For example the GKP vacuum is made out of S light-cone derivatives 
D + = D X {, i.e. we we can study this state by considering S excitations of type uf . Once the 
Bethe roots are found we can read of the anomalous dimension of the corresponding single 
trace through 



E = 2ig 2 (logQ 4 )' (i/2) + ig 4 (log Q 4 f (i/2) + c.c. + 0(g 6 



(74) 



where 



Qa(u) 



K a 

3=1 



U — U 



(«) 



)■ 



The algebraic equations are the Beisert-Staudacher equations 
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Figure 10: The several excitations around the BMN vacuum. 



Where the form of $(it) = u L + 0(g 2 ) and cr 2 (-u) = 1 + 0(g A ) is irrelevant for us while 



3=1 U J 



t <> =F i/2 



In (75) we colored blue the terms which can be set to one at one loop. In this case the Bethe 
equations take the usual universal rational form 



= n 

3 2"»1/ k,b^j,a u j "I. 



uf - + iM ah 



u[ a) -i5 n J ,fi„^ a) -4 6) -iM a6 



(77) 



where M a b is (one of the possible forms of) the PSU(2,2\A) Dynkin matrix [20]. 

Finally, let us end this section with a short useful dictionary between Lorents and £77(2) x 
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SU(2) indices, see e.g. [21]. We have 




D 2 




+2^ii 
+2^ii 





(78) 
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